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Abstract 

The equation for the conic sections describing the possible orbits in a potential V ~ r _1 is 
obtained by means of a vector constant of the motion differing from the traditional Laplace-Runge- 
Lenz vector. 
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The Laplace- Runge-Len z 1 ^^ 4 (or Hermann-Bernoulli-Laplace^) vector provides a simple 
and elegant way of obtaining the equation for the orbit of a particle moving under the 
influence of a radial inverse- square- law force.— This method also has the added advantage 
over the more common approach of direct integration of the equations of motion^ that the 
very existence of this single- valued constant of the motion explains the otherwise surprising 
degeneracy of the Coulomb /Kepler problem.— Our aim here is to show that there is another 
vector constant of the motion which makes a first principles derivation of the equation of 
the orbit even easier than the Laplace- Runge-Lenz vector or the direct integration methods. 

The force law F = —kr/r 2 (we assume k > 0; the repulsive case may be treated along 
the same lines) leads to the nonrelativistic equation of motion 



at r 



and to the well-known conservation laws for energy E = |mu 2 — k/r and angular momentum 
L = r x p. If we choose polar coordinates in the plane of the orbit, we have the usual 
expression L = mr 2 9, and also 



dp dv ■ dv L 

Tt =m Te Q = TeT 2 - (2) 



Combining Eq. (J2J) with Eq. (1) we obtain 

dv k „ 

Te=— L r - (3) 

Note that Eq. (3) implies immediately that v(0) = v(27r) and that the velocity vector traces 
out a circle if 9 is allowed to vary from to 2tt. Integration of Eq. (3) is trivial if we recall 
that dO/dO = -f: 

v-~0 = u, (4) 

where u is a constant vector. This vector is the constant of the motion we propose to use 
instead of the Laplace-Runge-Lenz vector in the derivation of the orbits. It is interesting to 
note that a quaternion equivalent to Eq. (4) was already known to Hamilton in 1845.— Un- 
fortunately, this simple result seems to have vanished from textbooks on classical mechanics 
after the first decade of the twentieth century.' 10 

The vector u may be evaluated in terms of the physical parameters of the problem. If 
we choose the angle so that the minimum value of r is at 9 = and adjust our coordinate 
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system so that v = t> y at r min , we have 

k 

u = («b (5) 

Because L = mr min t> , the energy is E = \mv\ — k/r min = \mv\ — kmvo/L. We solve this 
equation for vq and obtain 



k 

vo = t ± 



k\ 2 2E k 



+ — = -(l±e), (6) 



where e = Jl + 2EL 2 /mk 2 will turn out to be the eccentricity. The plus sign is required at 
r m ; n ; substituting this result into Eq. (5) provides us with an alternative form for u: 

u = ^ey. (7) 

Obtaining the equation of the orbit is now straightforward. By taking the scalar product 
of Eq. (4) with 6 and using v • — r6 = L/mr, y • = cos 9, we find 

L k k 

- = -ecos#. (8) 

mr L L 



The definition a = L 2 /mk allows us to write the solution for r as 



a 



(9) 



1 + e cos 9 

which is the usual equation of a conic section with one focus at the origin and eccentricity 
e. Note that the velocity v = /c[esin#? + (1 + ecos9)0}/ L = k[— sin#x+ (e + cos9)y]/L 
follows almost trivially from Eqs. (4) and (7). Hence the orbit in velocity space is always 
circular and characterized by v 2 + (v y — ke/L) 2 = k 2 /L 2 . More precisely, the orbit in velocity 
space is a circle of radius k/L and center at (0, ke/L) if the orbit in position space is a circle 
(e = 0), an ellipse (0 < e < 1), or a parabola (e = 1). On the other hand, if the spatial 
trajectory is a hyperbola (e > 1), the angle ranges from — cos~ 1 (— 1/e) to cos _1 ( — 1/e) 
only, and the velocity space orbit is a section of a circle in the upper v x ,v y plane with 
k(e 2 - l)/eL < v y < k(e + l)/L. 

The vector constant of the motion u is simpler (and its derivation certainly easie^) than 
the Laplace-Runge-Lenz vector, but the two are, of course, not independent constants of the 
motion. Indeed, the cross product of Eq. (4) with mL yields 

p x L — mkr = mu x L. (10) 



The left-hand side is the Laplace- Runge-Lenz vector A, so A = muxL. The three conserved 
vectors A, u, and L form a right-handed orthogonal system equivalent to x, y, and z after 
normalization. 

As a final comment of pedagogical interest, we present a second approach that may be 
useful in an elementary discussion of the Coulomb/Kepler problem. We begin with the 
expression for the energy E = \mv 2 — k/r. If we use v • = rO = L/mr to replace 1/r by 
mv • 0/ L in the potential term, we may write E as 

£ = L,,2_i.fl. (11) 
Z Lj 

If we complete the square in the velocity, we find 

,_, 1 9 mk 2 

B =2 mu -2F' (12) 
with u the same vector as in Eq. (4). It follows that the magnitude of u must be constant. 
To prove that this vector is conserved both in magnitude and direction, we must only take 
a derivative of u and substitute dv/dt = —kr/mr 2 and dO/dt = —Or into the result. The 
equation of the orbit may then be found as before. This approach does not require any 
integrations and amounts to a change of gauge from a scalar potential to an effective vector 
potential description of the problem. 
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